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SUMMARY 


Experience has shown that the flutter prediction process for airplanes can 
be greatly affected by strong concentrated nonlinearities which may be local- 
ized in the linking elements of the control mechanism, in the pivot joints of 
variable-sweep-wing systems, and in the connecting points between wing- and 
pylon-mounted external stores. The principle of equivalent linearization 
offers an efficient possibility for solving the related nonlinear flutter equa- 
tions in the frequency domain as a complement to the well-known time domain 
procedures. Taking as an example an airplane with nonlinear control character- 
istics, it is demonstrated how the equivalent linearization approach can be 
extended to rather complicated systems with multiple sets of strongly inter- 
acting, concentrated nonlinearities. 


INTRODUCTION 

Routine flutter analyses generally imply linearized representation of both 
the structural and the aerodynamic properties. This approximation has proved 
to be a useful basis for the flutter clearance of a large number of aircraft 
prototypes. There remains, nevertheless, a significant number of flutter cases 
suffering from rather poor agreement between analysis and test results. Many 
of these disagreements can be traced to structural nonlinearities. A survey of 
the various types of structural nonlinearities, their physical sources, and 
their effects on aircraft vibration and flutter is given in reference 1 , which 
indicates that strong concentrated nonlinearities are a common feature of the 
control systems of mechanically controlled airplanes. From reference 2, which 
presents a new experimental-numerical approach to determining the dynamic char- 
acteristics of hydraulic aircraft control actuators, it becomes obvious that 
flutter of aircraft with hydraulic controls may also be greatly affected by 
strong concentrated nonlinearities. References 3 and 4 focus on the special 
case of a modern variable-sweep-wing fighter airplane with concentrated non- 
linearities in the wing pivot mechanism and in the corresponding single-point 
external store suspension system. Reference 5 and, in particular, reference 6 
describe several concepts of how the governing equations of airplanes with con- 
trol system nonlinearities can conveniently be formulated in terms of consis- 
tent sets of both measured modal data and nonlinear force-deflection diagrams. 
The nonlinear flutter equations can be solved in the time domain by using ana- 
log computer techniques (see refs. 5 and 7) or by numerical integration. In 
addition to this time domain approach, promising attempts have been made to 
solve nonlinear flutter problems in the frequency domain by employing the prin- 
ciple of equivalent linearization (see ref. 8) . The effectiveness and accuracy 
of this equivalent linearization approach were impressively demonstrated for a 
semispan wing-aileron model with a single nonlinearity in the aileron hinge; 
the calculated and the wind-tunnel test results agreed very well (see ref. 1). 

Application of the equivalent linearization approach to systems with more 
than one nonlinearity creates some additional, though still solvable, difficul- 



ties. These difficulties are associated with an incompatibility between the 
input data representing the equivalent stiffness and damping properties of the 
nonlinearities involved and the corresponding output deflections, in coping 
with this problem, a recent investigation (ref. 9) describes the application 
of a method called the describing function method^ (ref. 10) to the special 
case of a flexible missile control surface with simple undamped free-play non- 
linearities in both the roll and pitch degree of freedom of the root support 
stiffness. 

The particular concern of the present study is the extension of the equiv- 
alent linearization concept to the flutter analysis of complete airplanes with 
strong hysteresis-type nonlinearities in the control system. Antisymmetr ical 
flutter of a sailplane involving strongly interacting rudder and aileron non- 
linearities is used as a realistic example to demonstrate the applicability of 
the method proposed. 


SYMBOLS 

A,B,C mass, damping, and stiffness matrices, respectively, defined in terms 
of physical deflections 

Aa,Ab,AC matrices of mass, damping, and stiffness changes, respectively, 
defined in terms of physical deflections 

b viscous damping coefficient of control surface hinge 

c stiffness coefficient of control surface hinge 

f frequency, 0)/2 tt 

F force or moment acting on control surface 

g absolute amplitude value, see equation (31) 

h bending deflection of quarter-chord line of lifting surface 

j imaginary unit, \f-T 

Z ha If -chord length 

Mf flight Mach number 

M,D,K generalized mass, damping, and stiffness matrices, respectively 

Am,Ad,Ak generalized matrices of mass, damping, and stiffness changes, 

respectively 


^A slightly modified form of the equivalent linearization approach. 
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N number of controls involved in flutter case 

P column matrix of external forces 

q column matrix of generalized coordinates 

Q column matrix of generalized forces 

R matrix of unsteady aerodynamic forces Rj. related to normal 

modes 

t time 

u column matrix of physical deflections 

V flight speed 

ot rotation about quarter-chord line of lifting surface 

3 control surface rotation about hinge line 

Y damping loss angle, 2C 

r diagonal matrix of damping loss angles 

e matching function, see equation (30) 

5 damping expressed as ratio to critical damping 

ri absolute amplitude value, see equation (29) 

A diagonal matrix of square values of circular normal frequencies 

Uj. = 2lTfp 

p air density 

T control surface chord length ratio (see fig. 11) 

(j) integration variable, (ot 

$ modal matrix of normal modes <I>r 

03 circular frequency 

Subscripts: 

A aileron properties 

L linear properties 

NL nonlinear properties 
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r normal mode index, r = 1, 2, . . n 

R rudder properties 

V,vi,a indices of concentrated nonlinearities involved 

Superscr ipts: 

P flutter speed 

T transposed matrix 

0 starting values 


NONLINEAR EQUATIONS OF MOTION 

Reference 6 offers a choice of several modal synthesis concepts which can 
conveniently be used to establish the aeroelastic equations of motion for the 
flutter analysis of airplanes with strong concentrated control System nonlin- 
earities. In accordance with one of these concepts (concept II of ref. 6) , the 
originally nonlinear airplane structure is physically converted to an artifi- 
cially linearized test configuration by replacing the nonlinear elements by 
linear stiffnesses with low damping. The normal mode characteristics of the 
linearized test configuration serve as a consistent basis for the calculation 
of both the unsteady aerodynamic reactions and a set of nonlinear coupling 
terms retransforming the test configuration to the actual system. The nonlin- 
earities can be determined statically in the form of force-deflection diagrams 
or dynamically by direct measurement of equivalent stiffness and damping values 
versus vibration amplitude. The equations of motion of the modified linearized 
test configuration, formulated in terms of physical deflections, can be written 
in matrix notation as follows; 

Au + B^u + Cj,u = P (1 ) 


where 

A mass matr ix 

Bl viscous damping matrix 

stiffness matrix 

P column matrix of external forces, for instance, unsteady aerodynamic 

forces 

u column matrix of physical deflections; u and u are first- and 

second-order differentials with respect to time t 

The dynamic behavior of the unchanged nonlinear system may be described by 

Au + Bu + Cu = P (2) 
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I 


where 


B - Bj, - ABl + 
C = Cl - ACl + 



(3) 


and where ACl and ABl denote the stiffness and damping properties of the 
artificial linear elements and ACjjl denote the amplitude-dependent 

stiffness and damping of the replaced nonlinearities. 


Development of the arbitrary deflection vector u in a series expansion 
of the normal modes 4>j- of the linearized test configuration yields 

u = $q (4) 


where 

$ modal matrix containing normal modes <I>j. as columns 

q column vector of generalized coordinates 

Substituting the modal transformation (eq. (4)) into equation (2), premultiply- 
ing by and taking into account equation (3) lead to the generalized equa- 

tions of motion of the unmodified nonlinear system 

Mq + (Dl - ADl + ADjjL)q + (Kl - AKl + ^KjjL)q = Q (5) 

where 

M = <I>Ta$ 

Dl = $'^Bl<I> 

ADjjl - 
Kl = 

Q = $Tp 

Consider the vth of N controls involved in a flutter case; the corre- 
sponding part of the modal matrix $ degenerates to the row matrix 

~ »Pv2'* * • / ^Vr ' * • *'^Vn^ 

with Bvr denoting the control rotation in the section where the control force 
is applied. Accordingly, the matrices ABj^l - ^Bl and degenerate 

to the 1 x 1 matrices 


- ABl)^> 


- AClI^* 


(6) 
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^L,V “ ^L,V = 

^NL,V “ ^L,V = CNL,v(M " °L,V 


( 8 ) 


where *3 l,v define the artificial hinge stiffness and damping of 

the vth control surface and Cnl^v^^v) ®nd bjjj,^ij(3v) define the amplitude- 
dependent stiffness and damping of the replaced nonlinearity of the vth control 
surface. Hence, 


N N 

(^NL,V " ^L,V) = '^v’’rt>NL,v(2v) “ t>L,v]^V 

V=1 V=1 


N N 

^NL " = ^T. (^NL,V “ ^L,v) = 2Z “^V^ , V < ^v) “ <=l,V^'^V 

V=1 V=1 


(9) 


The normal modes of the linearized test configuration satisfy the orthogo- 

nality condition, 

= M 


= AlM = Kl 


no) 


where 

M diagonal matrix of generalized masses 

. 2 
Kl diagonal matrix of generalized stiffnesses KL,r = “L,rMr 

diagonal matrix of square values of circular normal frequencies 

The generalized damping matrix Djj, which is not necessarily diagonal, was 
defined in equation (6) . Without damping coupling, matrix also becomes 

diagonal with the generalized damping elements DL,r* 

The unsteady aerodynamic forces P generally depend on time t, flight 
Mach number Mf, flight speed V, and air density p. Developing P in a 
series expansion of unsteady aerodynamic forces Rj. related to the normal 
modes 'J'j. leads to 

Q = R(Mf,V,P,t) q (11) 

where 


R(Mf,V,P,t) = [Rt,R 2,. . .,Rrr- • -»RJ (12) 

As mentioned previously, application of the equivalent linearization approach 
to nonlinear flutter problems requires a transformation of the differential 
equation (5) into the frequency domain. Accordingly, by assuming simple har- 
monic motions. 
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q(t) = 


(j = hv 


( 13 ) 


where u is the circular frequency, equation (5) reduces to 

[-(u2m + ju(Di, - ADl + ADul) + Kj, - AKl + A%l “ R (Mf , V,P ,u) ] q = 0 (14) 

Solutions of this equation can be obtained much more easily by expressing the 
viscous damping forces in terms of complex stiffnesses or damping loss angles. 
By so doing, equation (14) becomes 


’-(jo2m + j 


N N 

- > 'Yl,v ^Kl,v + 'V'nl,v(Pv) ^%L,V 

V=1 V=1 


+ Kl - AKl + AKnl - R(Mf,V,P,a))\q = 0 


(15) 


Damping can also be expressed by C as a ratio to the critical damping. The 
relation between C and Y is C = Y/2. In equation (15), T denotes the 
diagonal matrix of the damping loss angles Yr associated with the general- 
ized stiffnesses Kj-. The matrices Akj^ ^ and Akjjl v sre defined in equa- 
tion (9). The damping loss angles Yl,v^ Ynl,v(^v) represent the 

structural damping coefficients associated with the hinge stiffnesses Cjj^v 
and cjjL v(^v)f respectively, which are defined in equation (8). The matrices 
M, r, Ajj, and which describe the dynamic behavior of the modified linear- 
ized system, can be measured in a fairly simple ground vibration test (GVT) . 
These modal data and some related geometrical data are given in detail in 
appendix A for the sailplane taken as an example of a nonlinear system. 

Because of the high aspect ratio and the comparatively low maximum speed of 
this sailplane, incompressible strip theory is used to calculate the unsteady 
aerodynamic forces based on the measured mode shapes The method of deter- 
mining the nonlinear terms Ynl,v(^v) ^*^NL,v is described in the follow- 

ing section. 


EQUIVALENT LINEARIZATION APPROACH 

As is known from reference 8, an elastodynamic system with nonlinear 
stiffness and damping elements can be approximately described as a linear sys- 
tem for constant-amplitude vibrations at any arbitrary amplitude level. The 
fundamental idea of this equivalent linearization approach is based on the 
assumption that a nonlinear elastomechanical element can be approximately 
replaced by a linear substitute element with equivalent stiffness and damping 
energies when activated at equivalent amplitude levels. The accuracy of the 
approach, depending on the special problem to be investigated, can be assessed 
by procedures described, for instance, in references 11 and 12 for application 
to systems subjected to simple harmonic excitation. In addition, reference 12 
shows a simple way to solve problems with preloaded nonsymmetric nonlinearities, 
such as those arising in systems subjected to maneuver loads. 
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In accordance with reference 8 the equivalent linear coefficients of a 
nonlinear force-deflection diagram can be calculated from 


1 p2TT 

= — \ sin <j>) cos <1> d<|> 

1 

Ynl(P) r \ cos sin <j>) sin 4> d4> 

Cfjjj(6) 7 t 8 4j)=o 

where cnlCP) and Ynl(^) define the complex stiffness. 


(16) 


cJl(3) = Cnl( 3)[1 + j Ynl(2)1 (17) 

and where the force F is a nonlinear function of the deflection 6. Inte- 
gration is carried out over a full period of oscillation using (j> = wt as 
integration variable. 


In view of the particular flutter case to be dealt with subsequently, two 
special types of bilinear force-deflection diagrams, sketched in figures 1 
and 2, are evaluated by means of equation (16) . The diagram in figure 1 is 
characterized by a low stiffness ci for 

-6-, S 6 ^ (18) 

vdiere 6-] denotes the amplitude corresponding to the maximum stroke of the 
control surface. For 


-6-, > 3 > 3] (19) 

the stiffness assumes the much higher value C 2 because of kinematic limita- 
tion beyond the blocking point (see figs. 1 and 2). Hence it follows that 





. Ynl(^) = 0 


where 




arcsin 



(sin 2(|)-i 


+ 2<j>i) 





( 20 ) 


(21) 


Figure 2 illustrates a bilinear hysteresis-type force-deflection diagram. For 
amplitudes below the blocking point according to equation (18), the equivalent 
stiffness and damping values are 
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C2 

Cnl(^) = — - sin 2(j>i) 
2ir 


Ynl(^) = 


2 sin^ (j)-j 
2<j>-( - sin 2(|)-| 


( 22 ) 


where 


4>-| = arcc»s 1 - 


2Po 

C2Pi 


(23) 


where Fg is defined in figure 2. For amplitudes beyond the blocking point 
according to equation (19), the equivalent stiffness and dcunping values are 


C2 




°NL(^) = — (2'1>1 - 2(|)2 + it - sin 2 (J)i - sin 2<t>2) 
21T 


Ynl(3) = 


cos 24) I + cos 24>2 


24>i - 24)2 + ■'f ~ sin 24 >t - sin 24>2 


(24) 


where 


26 n - 3i 




4>t = arccos 


4)2 = arcsin — 


(25) 


where 3 q is defined in figure 2. 


SOLUTION OF THE NONLINEAR FLUTTER EQUATIONS 

Equation (15), which is usually written in the form of a complex eigen- 
value problem, can be solved for an arbitrary set of N equivalent stiffness 
and damping values °§l,v(^v) ^nd ySl,v(^v) which, because of their 
amplitude dependency, correspond to a definite set of deflections 3v^* Stan- 
dard flutter calculation techniques can be applied to determine the flutter 
boundary vrtiich is generally characterized by ein undamped harmonic oscillation 
of one of the generalized degrees of freedcm. The corresponding generalized 
eigenvector q^^ can be transformed into the physical deflections: 

viF = ^qr®" (26) 
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The deflections 3-^^ of the different control surfaces with nonlinear 
elements which are part of the flutter mode shape u^ can be determined by 

av®" = ^^v^r^ (V = 1, 2, . . N) (27) 

The deflections represent a consistent set of solutions if and only if 

the following condition is satisfied: 

= 0 (V = 1, 2, . . N) (28) 


where 


= l3v°l nv^ = |6v^l (29) 

In fulfilling condition (28) , the following iterative procedure may be 
employed. In a first step a set of amplitudes riv® corresponding to a set of 
equivalent stiffness and damping values and y8l respectively, 

is selected as input data for an initial flutter calculation. At the flutter 
speed resulting from this calculation, the related flutter mode deflections 

can be calculated by means of equation (27) . The difference between the 
initial amplitudes riv^ the flutter mode amplitudes both of which 

are defined by equation (29) , can be set into the more extended form 

N 

e = - Tlv°)^ = 0 (30) 

V=1 


where 


■ri Fm 0 

Mv ny 

9v^ = (V = 1, 2, . . ., N) (31) 

which means that the amplitude of the pith nonlinear element is kept constant 
for all further calculations; that is, 

(32) 

To determine an optimal change Ari'^ of the starting values riv® for the next 
iteration step, e is successively differentiated with respect to the ampli- 
tudes 

3e , 9ga^ 

= -2(gv^ - Hv°) + 2 - na°) - 

3riv° cr 3riv° 


(o — 1, 2, . . ., y— 1 , y+1 , . . . , N; 

V = 1, 2, . . ., y-1, y+1, . . ., N) (33) 
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The gradients 9e/3riv® can be approximately determined by successive small 
changes of the values ilv® resulting in N - 1 flutter calculations. To 
insure that each succeeding iteration step is going in an optimal direction, 
signs of the changes have to be chosen as follows: 


sign (Ari^) = -sign 



( 34 ) 


where the finite differences Ari^ are much larger than those used for the 
determination of the gradient 3e/3riv®. In approaching condition (30), it 
may be found that the difference - riv® changes sign. In this case, the 

next change Ati^, can be determined by interpolation between the two suc- 
cessive values of tIv^ ~ hv®- Consequently, the changes Ari^, for the next 
iterations have to be reduced to fractions of the last value of Ariv found 
by interpolation. The procedure as described above can be iteratively 
repeated until equation (28) is fulfilled (i.e., a matching point is 
obtained). Figures 3 and 4 give two examples of the functions riv^ “ hv® 
and e . 

Despite the comparatively large number of flutter calculations during the 
iteration process, the numerical effort in terms of computer time remains low 
because only one set of normal modes is required. Consequently, only one cor- 
responding set of unsteady aerodynamic forces is necessary. The applicability 
of the proposed approach is demonstrated in the following section for the exam- 
ple of a sailplane system with two complicated nonlinearities. 


APPLICATION TO CONTROL SYSTEM NCWLINEARITIES OF A SAILPLANE 
Description of Cases Studied 

To obtain a better insight into the mechanism of nonlinear flutter, three 
different configurations distinguished by different kinds of nonlinearities in 
the rudder and aileron control system of a sailplane are investigated. As a 
basis for selecting these configurations, some nonlinear data measured on the 
sailplane in dynamic and static tests are available. Thus, the nonlinear char- 
acteristics of the aileron system Ccin be derived from a statically measured 
force-deflection diagram (fig. 5) and from the dynamically measured aileron 
resonance frequency as a function of the amplitude 3^ (fig* 6) . The non- 
linear properties of the rudder system are available in the form of the 
statically measured force- deflection diagram in figure 7 with the standard 
trim stiffness removed. The special features of the three configurations can 
be described as follows: 

Configuration I .- To assess the importance of strong hysteretic damping on 
flutter behavior, damping is eliminated in this configuration. 
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Configuration II .- In this configuration the actual nonlinearities with 
measured hysteretic damping are taken into account, but the standard trim 
stiffness of the rudder system is eliminated. 

Configuration III .- The rudder and aileron system nonlinearities of this 
configuration are equivalent to those of the functionally complete sailplane 
including the rudder trim stiffness. 

The numerical values quantitatively describing the three configurations 
are detailed in appendix B. 


Results 

The nonlinear flutter boundaries of configurations I, II, and III are 
depicted in figures 8, 9, and 10 in the form of the amplitude ratios 
and Hr^/3t ,r as functions of the flight speed. Figure 8 shows that the arti- 
ficial configuration without damping results in a linear flutter behavior with 
a flutter speed independent of amplitude up to the blocking point at the kine- 
matic limit (riR^/3-| = 1). Above that amplitude, the flutter speed drops 

sharply as the amplitude further increases. This configuration has been inves- 
tigated to form a basis for demonstrating the effects of nonlinear hysteretic 
damping such as that inherent in the more realistic configurations II and III 
to be discussed next. 

As shown in figure 9, the influence of hysteretic damping results in a 
considerable stabilizing effect cxxnpared with configuration I. This means that 
between 195 and 235 km/hr, the nonlinear flutter boundary is characterized by 
increasing flutter speed as the amplitude ratios tIr^/3-1 ,r and iIa^/3-i 
decrease. However, the system at speeds within the above given range is stable 
only below a certain amplitude level, which can easily be exceeded by external 
excitation due to gust or maneuver loads resulting in violent divergent flutter. 
This special type of nonlinear flutter is contradictory to the wide-spread 
opinion that the effect of structural nonlinearities on flutter results always 
in limit cycle flutter vibrations. 

The functionally complete sailplane with a trim stiffness in the rudder 
control system is investigated in configuration III. As shown in figure 10, 
the additional spring stiffness results in considerable destabilization com- 
pared with configurations I and II. In a comparatively large speed range 
between 180 and 225 km/hr, external perturbations of not more than about 
50 percent of the blocking amplitudes of both the aileron and the rudder are 
sufficient to induce violent flutter. 

As for experimental verification of this flutter behavior, it is worth 
mentioning that seme time before the first routine flutter clearance on the 
basis of GVT and flutter calculations, a number of flight flutter tests were 
accomplished. The most remeurkable result of these tests was found at about 
180 km/hr. According to test pilot reports and observations from the ground, 
extremely lowly damped free vibrations with large amplitudes could be excited 
by pilot- induced rudder oscillations at an estimated frequency of between 3 and 
5 Hz, which forced the rudder to amplitudes near the blocking amplitude. On 
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the basis of figure 10, a small increment in speed would have resulted in 
violent flutter. It should be mentioned that the flutterlike vibrations were 
eliminated by mass-balancing the rudder. 


CCXfCLUDZNG REMARKS 

Based on standard flutter calculation techniques in the frequency domain 
and on the equivalent linearization approach, a method has been developed to 
predict the flutter behavior of complete airplanes with multiple sets of 
concentrated nonlinearities. The applicability of the method has been demon- 
strated for the example of antisymmetrical flutter of a sailplane with nonlin- 
earities in its control systems. The results are in good agreement with obser- 
vations during actual flight tests of the sailplane. 

In future investigations, emphasis should be placed on the following 
problems: 

1 . Investigation of service-life-dependent alterations of concentrated 

structural nonlinearities 

2. Investigation of mass, damping, and stiffness alterations due to pilot 

feedback 

3. Amendment of ground vibration test methods and flight and wind-tunnel 

flutter test techniques by paying more attention to nonlinear effects 

4. Development of suitable methods for the calculation of unsteady aerody- 

namic forces in the time domain for all flight speed ranges 

5. Application of digital and analog time domain techniques to solve tran- 

sient problems such as those due to gust loads, maneuver loads, and 
sudden failure of control system devices 

6. Investigation of the frequency-dependent dynamic properties of non- 

linear elements 

7. Investigation of nonlinear effects on design, test, and operation of 

flutter suppression and vibration reduction systems 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
February 8, 1980 
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APPENDIX A 


MODAL DATA OF A SAILPLANE IN A MODIFIED LINEARIZED TEST CONFIGURATION 

To obtain a largely linear test configuration, the nonlinear elements 
in the aileron, rudder, and elevator control mechanisms of a sailplane were 
replaced by linear lowly damped spring elements. A sketch of the sailplane 
investigated and the strip arrangement used to calculate the unsteady aerody- 
namic forces is shown in figure 11. Tables I, II, and III list the five lowest 
antisymmetrical normal modes and the geometrical data of the strip scheme. As 
shown in figure 11, the normal mode displacements referring to the midstrip 
sections are split up into the quarter-chord point bending deflection h, tor- 
sion a, and control surface rotation 6. The generalized masses Mj., the nor- 
mal frequencies £j. =0)j./2iT, and the damping loss angles Yr are listed in 
table IV. 

Since this investigation is of an antisymmetrical flutter case, only the 
rudder control system and the aileron control system must be taken into 
account. Correspondingly, in accordance with equation (7) , only the row 
matrices and describing the antisymmetrical hinge rotation angles 

^R,r ^A,r rudder and the aileron (V = r,a) are given in table V. 
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APPENDIX B 


QUANTITATIVE DETERMINATION OF NONLINEAR CHARACTERISTICS OF 
CONFIGURATIONS I, II, AND III 
Configuration I 

The following values are chosen for describing the nonlinearities of both 
the rudder and the aileron system by bilinear zero damping force-deflection 
diagrams such as those shovm in figure 1 : 

,A ~ N-in/rad c-] = 0 

C2^A ~ N-iu/rad <^2,R ~ 500 N-it(/rad 

3i ,A = 20° = 300 

The stiffnesses ,R» 02 ,r are estimated from figures 5 and 7, 

respectively. Stiffness c-j can be approximately calculated from 


M2 

Cl, A = — (2TA^)2(2TTfA,min)2 


(Bl) 


where 

M 2 generalized mass of normal mode r = 2 (see table IV) 

ratio of aileron chord length to wing chord length at inboard edge 
of aileron 

J half-chord length of the wing at inboard edge of aileron 

^A,min minimum resonance frequency of aileron system (see fig. 6) 

Application of equation (Bl ) implies that the normal mode r = 2 fundamentally 
consists of motion of the aileron system (see tables I, II, and III). The 
equivalent stiffnesses Crl^;^( 3) and Crl^r( 3) of the aileron and the rudder 
system calculated by means of equation (20) are shown in figure 12. 


Configuration II 

This configuration is characterized by hysteresis- type force-deflection 
diagrams in both the aileron and the rudder system. The force-deflection dia- 
gram of the rudder system shown in figure 7 can be approximated by a bilinear 
diagram of the kind sketched in figure 2. Because the trim stiffness is elimi- 
nated in this configuration, the following quantitative data can be selected: 
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APPENDIX B 


,R = 0 ,R = 30° 

°2,R “ 500 N-m/rad ®*o,R “ ^ N-m 

The corresponding equivalent stiffness and damping functions 
and Ynl,r( 2) plotted in figure 13. The equivalent aileron stiffness 

®NL,a( 3) can be calculated from 


M2 

®nl,a( 3) = — (2TA£)2[2rr fA(B)]2 (B2) 

The terms M 2 , and Jl are already defined with equation (Bl). Figure 6 

shows the resonance frequency £^(3) of the aileron system measured as a func- 
tion of 3» The equivalent aileron damping loss angle can be calculated 
approximately from 


and 


Ynl,a(B) 


4Fq,a 

■IT ®NL,a(^)^ 


YNL,a(2) 


4Fq,aPi ,A 
°nl,a<2)32 


(3S6i,a) (B3) 


(3 < 3i,a) (B4) 


vrfiere 


= 0.283 N— m 

The terms C(jl^A<3) and Ynl,a(3) are plotted in figure 14. 


Configuration III 

The aileron characteristics of this configuration are equivalent to those 
of configuration II. In contrast, the rudder characteristics change due to an 
additional trim stiffness attached to the rudder pedal mechanism. The equiva- 
lent stiffness can simply be calculated by adding the estimated trim stiffness 

Ct = 33.12 N-m/rad 

to the equivalent stiffness of configuration II; 


°NL,r(^) = ®NL,r(^) + °t (®5) 
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APPENDIX B 


For no damping in the trim stiffness, the damping loss angle changes to 




<=nl,r(P) 
°nl,r(^) + <=t 


Both c^^r( 3) and Ynl^r(®) 


are shown in figure 15. 


(B6) 


Artificial Linear Configuration 

Finally, the properties of the artificial linear elements were determined 
by means of dynamic tests to be 

''^LjA "■ 0.310 Yl,R " 0.272 

Cl^;^ = 51.25 N-m/rad Cl,R ~ 21.4 N-itv'rad 

Thus, all the data are known to set up the nonlinear equations of motion 
(eq. (15)). 
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TABLE I.- STRIP SCHEME AND FIVE ANT I SYMMETRICAL NORJIAL MODE SHAPES FOR WING 


strip no 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

Width, 

3, 

cm 

. 

. 

. 

. 

30 

157.5 

150 

103.5 

50 

50 

50 

50 

50 

59 

Half-chord 
of wing, 

length 
1 , cm 




46.38 

44.56 

41 .76 

39.44 

36.50 

33.50 

30.50 

27.50 

24.50 

21 .23 

Aileron 

ratio 

chord length 

.... 


0 

0 

0 

0 

0.264 

0.272 

0.280 

0.288 

0.296 

0 


h, 

cm 

. » 




0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Mode 
r = 1 

Otf 

rad 





0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


6, 

rad 





0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


h, 

cm 





-0.0094 

-0.0269 

-0.0263 

-0.0356 

-0.0363 

-0.0325 

-0.0250 

-0.0125 

-0.0060 

0.0060 

Mode 
r = 2 

Ct, 

rad 





0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


6, 

rad 





0 

0 

0 

0 

0.0513 

0.0520 

0.0528 

0.0535 

0.0560 

0 


h. 

cm 





0.0063 

0.0356 

0.0450 

0.0063 

-0.0375 

-0.0825 

-0.1340 

-0.1 950 

-0.2640 

-0.3470 

Mode 
r = 3 


rad 





0 

0 

0 

0 

0 

0 

0 

0 

0 

0 


6, 

rad 





0 

0 

0 

0 

0.0520 

0.0540 

0.0555 

0.0560 

0.0560 

0 


h. 

cm 





0.0281 

0.1710 

0.3330 

0.2720 

0.1250 

-0.0188 

-0.1980 

-0.4420 

-0.7360 

-1.110 

Mode 
r = 4 

a, 

rad 





0 

0 

0 

0 

0 

■ 

■ 





B, 

rad 





0 

0 

0 









h. 

cm 





0 

0 


0 


H 





Mode 
r = 5 

a. 

rad 





0 

0 


0 


m 



H 



B, 

rad 


• 



0 

0 


0 


■ 



IB 
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TABLE II.- STRIP SCHEME AND FIVE ANTISYMMETRIC NORMAL MDDE SHAPES FOR ELEVATOR 


Strip no 

1 

mgm 

3 

4 

5 

Width, 

s, cm 

mm 


30 

30 

30 

Half-chord length 

of elevator, S-, cm 

25.75 

24.5 

23.25 

20.81 

19.69 

Mode 
r = 1 

h, cm 

0 

0 

0 

0 

ililiR 

a, rad 

0 

0 

0 

0 

0 


h, cm 

-0.0125 

1 

o 

• 

o 

00 

CD 

-0.0513 

-0.0806 

-0.0844 

■an 

a, rad 

0.0008 

0.0008 

0.0013 

0.0017 

0.0017 

— 1 


0.0575 

0.1725 

0.3581 

0.5406 

0.7125 


a, rad 

-0.0005 

-0.0009 

-0.0005 

-0.0002 

-0.0009 


h, cm 

-0.0125 

-0.0238 

-0.0531 

0.0738 

-0.0863 


a, rad 

0 

0 

0 

0 

0 


h, cm 

0.0838 

0.2495 

0.4969 

0.7431 

1 .0000 

H 

a, rad 

0 

-0.0005 

-0.0012 

-0.0013 

0 









































TABLE III.- STRIP SC3iEME AND FIVE ANTISYMMETRIC NORMAL MODE SHAPES 


FOR VERTICAL TAIL 


Strip no 

1 

2 

3 

4 

5 

Width, 


50.64 

44.73 

37.00 

28.60 

22.10 

Half -chord length of 

vertical tail, il, cm . . . . 

20 

30 

40 

31 .5 

21 

Rudder chord length 

ratio, 


0.450 

0.459 

0.469 

0 



0.0313 

-0.0940 

-0.1565 

-0.0940 

-1.1200 

Mode 
r = 1 

ct, rad 

0.0043 

-0.0043 

-0.0113 

-0.0140 

0.1020 


3, rad 

0.1000 

0.1000 

0.1000 

0.1000 

0 


h, cm 

-0.0625 

-0.0838 

-0.1163 

-0.1525 

-0.1838 


a, rad 

0 

0 

0 

0 

0 

H 

3, rad 

0 

0 

0 

0 

0 



0.2688 

0.3975 

0.6063 

0.8563 

1.044 

Mode 
r = 3 

rad 

0.0025 

0.0025 

0.0032 

0.0061 

0.0085 


3, rad 

0 

0 

0 

0 

0 


h, cm 

-0.2038 

-0.2288 

-0.2750 

0.3400 

-0.3688 

Mode 
r = 4 

otf rad 

-0.0018 

-0.0018 

-0.0025 

-0.0039 

-0.0025 


3, rad 

0.0024 

0.0024 

0.0024 

0.0040 

0 


h, cm ........ 

-0.2063 

-0.1325 

-0.0250 

0.0938 

0.1212 

Mode 
r = 5 

0(, rad 

-0.0017 

-0.0019 

-0.0019 

-0.0012 

0.0058 


3, rad 

0.0064 

0.0064 



0.0064 

0.0064 

0 
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TABLE IV.- MDDAL CHARACTERISTICS 

NORMAL MODES 


OF THE FIVE LOWEST ANT I SYMMETRICAL 
OF A SAILPLANE 


Mode , r 

Mj. , kg-cm^ 

fr. Hz 


Remarks 

1 

19.22 

1 .68 

0.272 

Rudder 

2 

3.50 

4.90 

.310 

Aileron 

3 

12.90 

5.39 

.042 

Tail - aileron 

4 

15.40 

6.69 

.019 

First antisymmetrical wing 
bending 

5 

3.24 

L 

8.57 

.078 

Elevator 


TABLE V.- ROW MATRICES CONTAINING HINGE ROTATIONS 6R,r ^ND OF RUDDER 

AND AILERON ACCORDING TO EQUATION (7) 

= To. 1 ,0,0,0.0024,0.0064] 

= [0,0.0573,0.0520,0,0] 
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V, km/hr 



Figure 3.- Some characteristic functions near matching point. 
Configuration II; = 0.858. 



0 


V , km/hr 



Figure 4.- Functions ^ versus tia® special case of 

three different matching points. Configuration III; riR®/3i n = 0.5. 



Aileron 
hinge 
moment, 

F(3^) , N 



Aileron hinge angle. 


6^, deg 


Figure 5.- Measured force-deflection diagram of aileron system for two different 

amplitudes. Antisymmetr ical case. 




Aileron 

resonance 

frequency, 

f^(3) , Hz 



Figure 6.- Measured resonance frequency of antisynunetr ical 
aileron vibration versus hinge angle. 










Figure 9.- Nonlinear flutter boundary for configuration II. 
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Rudder actuating force 



Figure 17.- Schematic view of a sailplane. 






Figure T3.- Hinge stiffness and damping of rudder without trim stiffness 
versus hinge angle. Configuration II. 




Figure 14.- Aileron hinge stiffness and damping versus hinge angle 

Configurations II and III. 
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Damping 
loss angle, 


Y ( 3 ) 


angle, 3 , deg 

K, 


iping of rudder with trim stiffness 
Configuration III. 
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